Reading charge transport from spin dynamics on the surface of topological insulator 
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Resolving the conductance of the topological surface states (TSSs) from the bulk contribution has 
been a great challenge for studying the transport property of topological insulators. By developing a 
non-purturbative diffusion equation that describes fully the spin-charge dynamics in the strong spin- 
orbit coupling regime, we present a proposal to read the charge transport information of TSSs from 
its spin dynamics which can be isolated from the bulk contribution by time-resolved second harmonic 
generation pump-probe measurement. We demonstrate the qualitatively different Dyaknov-Perel 
spin relaxation behavior between the TSSs and the two-dimensional spin-orbit coupling electron gas. 
The decay time of both in-plane and out-of-plane spin polarization is naturally proved to be identical 
to the charge transport time. The out-of-plane spin dynamics is shown to be in the experimentally 
reachable regime of the femtosecond pump probe spectroscopy and thereby we suggest experiments 
to detect the charge transport property of the TSSs from their unique spin dynamics. 
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PACS numbers: 73. 25.-1-1, 72.25.Rb, 72.10.-d 



Topological insulators (TIs) are a class of time-reversal 
invariant matteri^. Recently, both theorj*^"— and 
experiments^^— focus on the transport property of the 
TIs which have spin-momentum locking surface states. 
The helical topological surface states (TSSs) break the 
Fermi doubling and have very strong spin-orbit coupling 
(SOC) which gives their unique charge and spin transport 
propertiesi^ii^. However, in most Tis the bulk states also 
contribute the conductance, which makes it a great chal- 
lenge to extract the transport properties of the TSSs, 
such as mobility, from the experimental data^i^i. On the 
other hand; the spin polarization of the TSSs has been 
detected by the time-resolved second harmonic genera- 
tion pump-probe measurement which can separate the 
surface response from the bulk^i^. This has not been 
possible until now because the complete spin dynamical 
theory of TSS is still absent due to the strong SOC of 
TSSs which can not be treated perturbatively as it is 
usually done in the traditional diffusive equation. 

Here, we developed a non-perturbative spin dynamic 
theory to fully capture the spin dynamics of the TSSs 
which reveal the qualitative difference to that in other 
SOC 2DEGs such as GaAs quantum welli^^il. The spin 
polarization along x,y,z direction are decoupled to each 
other even though the TSSs have strong SOC. The decay 
time of both in-plane and out-of-plane spin polarization 
is predicted to be exactly equal to the charge transport 
time, Ttr, instead of being proportional to the inverse of 
rig^Tp, which is the commonly label Dyakonov-Perel (DP) 
mechanism. Here is the SOC frequency and Tp is the 
momentum scattering time. We estimate the time evo- 
lution of the out-of-plane spin dynamics which we show 
to be in the experimentally reachable regime of the fem- 
tosecond pump probe spectroscopy^. Because the surface 
spin dynamics can be isolated from the bulk contribution 
by recent second harmonic pump probe measuremenli^i^, 
we propose to read charge transport property of TSSs 



from the spin by the pump probe measurement. 

Non-purturbative spin-charge dynamic equation on the 
surface of TI - We assume that the low energy TSS has 
the Dirac-like Hamiltonian plus nonmagnetic short range 
disorder potcntialii^ 



(1) 



where v is the constant velocity, z is the unit vector per- 
pendicular to the surface of TI, k is the wave vector of 
the electron, and (Tq is 2 x 2 identical matrix. Here we 
assume ft = 1 and consider only spin independent scat- 
tering. The velocity operator has the form 
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dH 
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and is proportional to the in-plane spin polarization per- 
pendicular to the velocity direction. 

The dynamics of the spin-charge polarization, as a non- 
equilibrium process, can be described by quantum Boltz- 
mann equatior 
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where Ho{k) = v{z x cr) ■ k and 
g{ek,R,t) = 



(3) 



(4) 



is the angular distribution functionSS^ is the 2x2 ma- 
trix of spin-i Kcldysh Green's function, fc/ is the Fermi 
wave length. 
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is the density matrix, Tp is the momentum scattering time 
at Fermi surface, N{€k) is the density of states, rii is 
the 2D density of the nonmagnetic impurities and 9 is 
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FIG. 1. (a) Spin polarization configure at ky — 

cliannel in 2DEG. The solid(hollow) circle means creat- 
ing(annihilating) an electron into the system which is cor- 
responding to C^(C) operator. The clockwise SHM is con- 
structed by 'i+q t''*.i^"''I ^Cfe+5,-t-(the red creation and an- 
nihilation operators) and the counter-clockwise SHM is con- 
structed by c^_q^^Cfc4 + cl.j_Ck-q, tithe blue creation and 
annihilation operators). Here t i-l) as the subscripts indicates 
spin along -|-(— ) y direction. The collective spin polarization 
plot in the rectangle boxes is in a: — jz plane, (b) Spin po- 
larization configure at fcy = channel on the surface of TI. 
Because there is only one spin sub-band on the Fermi surface, 
the collective spin polarization arises from the summation of 
the half clockwise SHM and half counter-clockwise. 



the angel between kf and the axis. Eg. IA9I contains 
both commutation and anti-commutation relation of g, 
H , etc. and can only be solved perturbatively. However 
the difficulty of obtaining the spin dynamics equation 
in the strong SOC regime lies on the fact that the spin 
precession angle 2vkfTp, is too large to be treated per- 
turbatively. This fact makes it difficult to apply Eq. IA9I 
directly to the strong SOC regime. To circumvent this 
problem, we multiply Cj where j = 0, x,y, z on both sides 
of Eq. IA9I and calculate the trace. Using the fact that 
Tr((TjfTk)/2 = Jjk, the quantum Boltzmann equation can 
be written in the classical spin-charge 4D space after in- 
tegrating out E asiS 



(6) 



where gj{pj) ~ Tr[g(/5)(Tj/2] (the explicit form of K is 
shown in the supplementary material). Here, we aban- 
don the idea of gradient expansion in terms of VlTp and 



derive the spin-charge diffusion equation by simply multi- 
plying on both sides of Eq. (|Alip and integrating out 
the angle 9. To simplify our discussion, we assume the 
uniform spin polarization along y direction, say dyp = 0. 
In this case, using the relation Ea. (|A8[) . we obtain the 
non-purturbative spin-charge dynamic equation 



Pj = DjkPk, 



(7) 



where D — J ^K^^ and K^^ has a complicated but 
analytically and numerically tractable form shown in the 
online supplementary material. 

Eq. (jA13[) is a generalized spin diffusion equation and 
valid from the weak to the strong SOC regime. For the 
spin-momentum locked states, we have found that the 
spin polarizations along x, y and z direction are com- 
pletely decoupled even though the TSSs have very strong 
SOC. As a result, spin helical modes (SHMs) are ab- 
sent on the surface of TI. The absence of SHM is due 
to the fact that for the TSSs there is only one spin de- 
gree of freedom on the Fermi surface. For an intuitive 
understanding of this fact, in Fig. [TJ we schematically 
show that at the ky = channel, the spin polarization 
modes of both 2DEG and the TSSs. For the spin po- 
larization of the 2DEG shown in Fig. [TJa), there are 
two spin sub-bands on the Fermi surface which allow 
two SHMs with different two relaxation time a^^'^^i^^ . In 
Fig.[TJb), because the non-equilibrium state is around the 
Fermi surface and Ef 3> fcT, the spin sub-band below the 
Dirac point is still completely filled. Accordingly, at the 
ky ^ Q channel the terms such as X]/c>o ^fc i^k+q.^ and 

Sfe<o t^'^'-'-' ^^^'^h ^^'^ electrons to the spin sub- 
band below the Dirac point and thereby will not con- 
tribute to spin polarization. Here t {\.) as the subscript 
of the creation or annihilation operator indicates spin 
along +(— ) y direction. Therefore, the non-equilibrium 
spin polarization of TSSs around Fermi surface has only 
a z-component which is built from half clockwise terms 
containing '^k>o^k+q-\^k.i and half counter-clockwise 



The detailed discus- 



terms containing X]fc<o ^\ 
sion is presented in the online supplementary material. 

In plane spin dynamics along x direction.- The spin 
polarization is completely decoupled for the spin polar- 
ization along x-direction. Its non-purturbative dynamic 
equation determining its complex frequency takes the 
form 
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=0, 
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where I is the mean free path, 0=1- 
and ^x{y) = l<lx{y)- Here we have Fourier transformed dt 
and dx(y) to —iui and iqx{y) which are the frequency and 
wavelength of the spin polarized wave. The in plane spin 
polarization of the TSSs has a very unique property: it is 
proportional to the charge current perpendicular to it on 
the TI surface. Therefore the charge current relaxation 



time must be equal to the in plane spin relaxation time. 
A valid spin dynamic equation of the TSSs should be 
able to reflect this fact. To keep this argument in mind, 
let us focus on the spin polarization along x-direction. 
First we assume that <ti ^so and A^; <^ i^so- The 
eigen-frequency of the spin dynamics takes the form 



1 
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which can be Fourier transform to the real space as 



dtPx = -Dsdlp^ 



2r„ 



(10) 



where = WjTp/2. The prior theory of spin dynam- 
ics of the TSSs (Refi) is based on the Kubo formula 
and needs to expand the spin-charge kinetic equation in 
terms of uTp and A^; to their leading order. Making the 
same approximation, our dynamic equation of the spin 
polarization along x direction becomes 



A?. 

iujTp = 1 + = 0, 



which is equivalent to 



OtPx = --Dsd^Pa: . 

Z Tp 



(11) 



(12) 



Eq- (|12X is consistent with the spin diffusion equation 
in Ref. i3| obtained from the linear response theory in 
the weak SOC limit. However we have found that when 
w = 0, Eq. (HI]) gives A^^ = g^P = -4 which is of the 
same order but four times larger than the exact solution 
A^ = — 1. Both Ax = —4 or A^. = — 1 are inconsistent 
with the assumption that \Ax\ <C 1. On the other hand, 
the non-purtabativc result of Eq. [10] indicates the spin 
relaxation time is Ts = 2Tp, instead of Ts = Tp indicated 
by Eq. [12] based on the purtabative calculation. It is well 
known that charge transport time Ttr for steady current 
in the system with Dirac-like Hamiltonian is double of the 
momentum scattering time Tp due to the spin-momentum 
locked states^. Therefore, we our non-perturbativc the- 
ory systemically proves that that the spin relaxation time 
and charge transport time rtr are equal for the TSSs, 
Ts = Ttr. This conclusion is consistent with one of the 
unique properties of the Dirac-like Hamiltonian, that the 
spin operator is proportional to the velocity operator. 

In the case of uniform in plane spin polarization, with- 
out assuming fi <^ ^so, there is another solution of 
Eq.\M6\ 



(13) 



which gives a damped oscillatory spin dynamic mode. 
Due to the strong SOC, this spin oscillation is accompa- 
nied by an AC current with the same frequency. This AC 
current may provide another way to detect the transport 
property of TSSs. 

Out of plane spin dynamics.- next we focus on the spin 
polarization perpendicular to the surface. The dynamic 
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FIG. 2. (a) The second harmonic pumping probe measure- 
ment is proposed to observe the dynamics of the out-of-plane 
spin polarization. The circularized pump light with the fre- 
quency cJo is followed by the linear polarized probe light with 
the same frequency. The time delay is t. The linear polarized 
reflected light will only contain the information of the surface 
states^. 6 is the rotation angle of the linear polarization due 
to the Kerr effect, (b) The evolution of the uniform out of 
plane spin polarization. According to our theory, take the ex- 
perimental parameters from the sample Q2 in table 1 of the 
Ref.i^ as kf = 0.032A"\ Ej = 84meV and kfl = 69. 



equation of the spin polarization along z direction takes 
the form From Eq. (|A16p 



1 



f22 + Al 



= 0. 



(14) 



When ujz = for the steady spin polarization, we have 



Ax = i\ 2 + n 



2kfl, 



(15) 



which means the diffusive length of out of plane spin is of 
the same order as the Fermi wave length. This is in con- 
trast to the diffusive length in other materials such as the 
GaAs/AlGaAs where it is of the same order of (or much 
smaller) then the mean free path I in the strong (weak) 
SOC regimei^iiii^^i^. When the spin is uniformly polar- 
ized perpendicular to the surface, the eigenfrequency has 
the form 



1 -4f72 



1 



(16) 



where the approximation is valid when Qso S> 1. Eq. 1161 
indicates that the spin along the ^-direction will os- 
cillate in a damped form with the frequency equal to 
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riso = 2vkfTp and the relaxation time given as Tg = rtr, 
which is the same to the decay rate of the in plane spin 
polarization. This provides a reliable way to measure 
the charge transport time Ttr- Because the spin along 
z-direction is not coupled to any other spin and charge 
density, this oscillation will not induce current oscillation 
or charge density oscillation. The experiments have been 
able to freely tunc the Fermi surface of TI close to or away 
from the Dirac point in a wide range^^. A circular polar- 
ized pump beam can generate spin non-equilibrium polar- 
ization of TSSs which has been predicted theoretically^^ 
and verified experimentally^!. The spin dynamics of the 
TISS can be detected by the time-resolved second har- 
monic optical pump-probe measurements. As only TSSs 
contribute the second harmonic generatioit^, this type 
of measurement can isolate the surface response from the 
bulk. If the Fermi surface is turned to 84 meV— above 
the Dirac point, the experiment should observe about 24 
THz oscillation through the Kerr rotation angle. The 
femtosecond pump probe spectroscopy can be used to 
measure this oscillation. At the same time, the decay 
rate of this spin oscillation is Ttr which gives us the charge 
transport time of the TSSs (FigH]). 

In conclusion, we develop a non-purtabativc method to 



obtain the generalized spin dynamic equation of the TSSs 
which is in the strong SOC regime. Based on this equa- 
tion, we have found the in-plane spin relaxation time of 
the pure exponential decay mode is exactly equal to the 
charge transport time which is consistent to the unique 
property of the Dirac-likc spectrum considered, i.e., the 
charge current is proportional to the in plane spin polar- 
ization. This is a key check of the validity of our theory. 
We also predict two fast oscillatory modes of spin po- 
larization perpendicular and parallel to the surface of a 
topological insulator, which cannot be obtained from the 
prior spin diffusion equations, based on perturbative ap- 
proaches. Therefore the generalized spin diffusion equa- 
tion is necessary to quantitatively understand the spin 
dynamics. At last, wc have shown how to read the charge 
transport properties of TSSs from the out of plane spin 
dynamics which may provide an alternative way to de- 
tect the conductance of the TSSs isolated from the bulk 
contribution. 
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SUPPLEMENTARY ON-LINE MATERIAL FOR "READING CHARGE CURRENT FROM SPIN ON 

THE SURFACE OF TOPOLOGICAL INSULATOR" 

In this appendix, we present the detailed calculation of the spin dynamics on the surface of topological insulator(TI) 
and how the spin polarization is constructed by the non-equilibrium electrons. 



Appendix A: Spin-charge dynamic equation 

We assume that the low energy topological surface state (TSS) has the Dirac-like Hamiltonian plus nonmagnetic 
short range disorder potential^>^ 

H = v{zxcT)--k + VaJ:,6{r-r,)ao, (Al) 

where v is the constant velocity, z is the unit vector perpendicular to the surface of TI, k is the wave vector of the 
electron, and ctq is 2 x 2 identical matrix. Here we assume H = 1 and consider only the spin independent scattering. 
The velocity operator of this type of Hamiltonian has the form 

and is proportional to the in plane spin polarization perpendicular to the velocity direction. 

The dynamics of the spin-charge polarization, as a non-equilibrium process, can be naturally described by quantum 
kinetic equationSS 

dtG'' + Vr • [i{Vfc, G^}] + *[i/o(fc), G^] 

(S^G^ _ G^E^) - (G^E^' - E^G^)j , (A3) 

where t and R are time and real space variants, is the 2x2 matrix of spin-^ Keldysh Green's function, V = v{zx(t) 
is the velocity operator, Ho{k) = v{z x cr) ■ k and Yl^^"^'^^ is the retarded (advanced, Keldysh) self energy due to 
disorder potential. The term [5{Vfc,G^}] = J is natural definition of current matrix. The third term in the Eq. 
(|A3[) is the spin torque exerted by the Rashba SOI. In this work, the operator with a hat means it is a matrix. In the 
equilibrium state, the Keldysh Green's function satisfies^ 

G„^ = (G«_G^)tanh(|--^), (A4) 

where Gq stands Keldysh Green's function in equilibrium, fc^ is the Boltzmann constant and T here is the system 
temperature. Based on Eq. (jA4p and we assume that the Keldysh Green's function generally has the form 

G^^' = -2m5{E - e+)h+ - 2Tii5{E - )/i", (A5) 

where — ±vk, h^{k, R^T) is the distribution function, defined as 

/°° HP 
^^Gl^iR^t) (A6) 

and ± label the distribution function of the electrons spin sub-band abovc(bclow) Dirac point. Normally the non- 
equilibrium distribution function 5h locate around Fermi energy. Therefore we define the non-equilibrium thermal 
average distribution function as 

g^{ek,R,t)= j N{ek)dek5h^{k,R,T) (A7) 

where = k/k is the unit vector along fc direction and N{e^.) = k/hv is the density of state. In this work, we 
consider the case that the Fermi energy is above the Dirac point and much larger than fc^T. 

To obtain the kinetic equation of distribution function, we introduce the non-equilibrium density matrissi^ 

i f (Pk - f d9 
= 2^^ J J2^'' - J ^^^(%'I^'t)' 
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where 9 is the angel between kf and x axis. By integrating out energy E and The Eg. IA3l can be converted to 

^tg + WR■{\vr9} + ^mkf)r9] + - = ^-^^ . (A9) 

where 



is the momentum scattering time at Fermi surface and rii is the 2D density of the nonmagnetic impurities. 
The non-equilibrium thermal average distribution function and density matrix can be generally written as 

9 = 9c(^l3 + 9x(^x + 9y<^y + 9z(^z, 



(AlO) 



If we multiply where i = 0,x,y, z on both sides of Eg. I A3 1 and calculate the trace, using the fact that Tr((Ti(7j )/2 — dij 
the quantum kinetic equation can be written in the classical spin-charge 4D space after integrating out E as^ 



K 



( 90 

gx 

9y 

\ gz 




(All) 



where 



K 



Q 


-iAy 


iA 


-iky 


n 










n 



\ r2soCos(0) l^so sin(0) 



- \ 

-flso cos{9) 

-n^o sm{9) 

n ] 



(A12) 



r2 = 1 — iwTp, fiso = 2vkfTp and A^-jj,) — lqx(y) is from the anti-commutator on the left hand side of Eq. (|A9P 
which gives the spin-charge coupling. Here we have Fourier transformed dt and to —iuj and iqx(y) which are 

the frequency and wavelength of the spin polarized wave. For convenience, in the rest of this paper we will express 
frequency in units of l/r^ and wave vector in units of 1/Z, unless otherwise stated. The difficulty of obtaining the spin 
dynamics equation in the strong SOC regime lies on the fact that the three terms, spin precession angle 2vkfTp, the 
dimensionless time and spacial spin relaxation rate ujTp and g/, may be large simultaneously and cannot be treated 
perturbatively. Here, we abandon the idea of gradient expansion in terms of uiTp and ql and derive the spin-charge 
diffusion equation by simply multiplying on both sides of Eq. (jAlip and integrating out the angle 9. To simplify 
our discussion, we consider the spin wave vector is only along x direction and takes qy = 0. In this case, using the 
relation Eq. (jAll[) the spin-charge dynamic equation of the density coefficient Pc(x,y,z) can be obtained as 



9a \ 




I Pa 


9x 




Px 


9y 






9z 1 




\ Pz 



d9 
2^ 



where D ~ f ^ and K ^ takes the form 

/ + n^^^ i cos{e) sm{9)Axnl^ 

i cos{e) sin(6i) A^fi2^ Cl (jl^ + Al + sm^{9)nl 
-iAx (n^ + cos2(6i)f22^) -f2cos(6i) sm{9)n^,^ 




(A13) 



V 



id sm{9)Ax^s 



cos{9) {n^ + Al)n 



'iAx [n^ + cos^i9)n'^, 

-ncos{9) sm{9)n^^^ 

+ cos^{9)n^^n 



—iClsm{9)Axd=^ \ 

cos{9) (n^ + Al) n 
sm{9)ns. 
n[n^ + Alj J 



f!2 + 02jf!2+A 



^2 (n^ + cos^{9)ni 



(AM) 



The Eq. (|A13p is a generalized spin diffusion equation and valid from the weak to the strong SOC regime. It is noted 
that the denominator of K^^ is the function of cos2(0). By using the property of trigonometric functions without 
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calculating the integral of 0, it can be easily proved that the spin polarization along x and z direction is not coupled 
to other components. Therefore the Eq. (|A13p can be block diagonalized as 



1 



A, 



l-f^lxj-^ ^ r- = 

'^i Aiy/(V+Aiy/n-^ +Ai+ni 



Po 
Py 



= 



(A15) 



Al^/sP+nl^ 




A2 







+nj^ y/iV +Ai+ni^ 



(A16) 



For the spin-charged coupled equation, Eq. IA151 when the spin is uniformly polarized, say A^; = 0, the off diagonal 
terms in the spin dynamic matrix vanish. Therefore, the charge and spin densities for the uniform polarization are 
independent. 



Appendix B: In-plane spin dynamics under diffusive approximation 



Now, let us focus on the in-plane spin polarization along x direction. It has been proved exactly to not couple the 
spin polarization in other directions. Therefore we only need to calculate the element D22 which takes the form 



D22 = 



d0 



- sin^ 9 cos eK^h"^^ + + AI + sin^{9)n 



2tt 



(Bl) 



Due to the fact that sin^(0 + tt) cos(0 + n) ~ — sin^ 9 cos 9, we can drop the fist term in the integrated function of 
Eq. IBll and have 



d9 



n ( + Al + sm^ {9)nl 



2tt 



02 + nj) 02 + A2 (172 + cos2(6')f22 



(B2) 



Taking the diffusive approximation <C 1 and uiTp <C 1 while flso 3> 1, we expand the Eq. IB2I up to the first 
nonzero order in terms of A^. and f2,n as 



d9 



n m2 + A2 +sin2(6i)f72 



+ fi2^j fi2 + A2 (02 + C0s2(6')fi2^ 

d9 nsm^9 1, A2 

-= = ^ -(l+tLOT -). 

27r 1]2 + A2 cos2 61 2^ 2 ^ 



Substituting Eq. IB3I to the generalized spin dynamic equation Eq. IA13[ we have 

1 ~2 

lUJTp = 1 + -A^, 

which can be Fourier transformed to the real space as 

dtPx = -Dsdlpx - 



(B3) 



(B4) 



(B5) 



where Dg = v^Tp/2. This is equivalent to the in-plane spin diffusive equation in Ref|^. (Here we assuming spin wave 
vector along x direction and the term contain dy is zero). 
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Appendix C: The spin polarization texture in 2DEG and on the surface of TI 

111 this section, we intuitively show why spin helix mode is absent on the surface of TI. In the following derivation, 
we focus on the ky = channel to simple our calculation. 
The spin wave operator in real space is defined as 



and 



First, let us consider the spin polarization on the 2DEG. The operator 5*+ actin on electron states gets an spin 
density wave which has the form 



= Yj V^fc(l - hk) {cos{qx)(T^ - s\n{qx)ax) ■ (C3) 



where /ifc,t(4.) is occupied probability which is the diagonal element of the distribution function in Eq. IA6I Here, in 
the case of g ^ fc/, we make an approximation that hk+q.f = hk.i- The spin polarization given by Eq. IC3I in three 
directions are (Sz) = cos{qx), (Sx) = — sm{qx) and {Sy)^0 which is a contour-clock wise spin helix mode in a; — z 
plane. 

Similar the clockwise spin helix mode is constructed by 

5_ = {S.{x)) = J2 (v/A-.,tA,i(i - A-«,t)(i - fk,i)^ ( y Zl 

e-*9^ / 1 i 



i -1 



« E V^fc(l - hk) (cos(gx)cr^ + sm{qx)ax) , (C4) 

A: 

whose spin polarization in three directions are (Sz) = cos{qx), (Sx) ~ siii(ga:) and (Sy) ~ 0. 

On the surface of TI, there is only one spin-sub band at the Fermi surface. The other sub-band is below the Dirac 
point and thereby far below the Fermi energy in the condition Ef ^ kT. Accordingly, for the counter-clockwise 
spin helix mode in Eq. IC3[ hk+q,-\ and /ife,4, can not be approximately equal because they belong to different spin-sub 
band and has the energy difference about 2Ef. For example, when fc < 0, as shown in Fig. lb in the paper, the state 
\k + (jf,t) is below the Dirac point and completely occupied so that ft,fc+q<o,4. = 1. Similar we also have hk>o.f = 1. 
As a result, the terms containing J2k>o ^1 i^*:+9,t i^k+q<o ^l+q ^^k,i) becomes zero because they are proportional 
to ^/l — ^fe+(j<o,t = 0(-\/l — hkyo.i = 0). This implies hat only half terms in Eq. IC3I which have the form 

k + q>0 k<0 

A;+g>0 ^ ^ fc<0 ^ ^ 

will contribute to the spin polarization. Similar, for the clockwise spin helix mode, only half terms in Eq. IC4I which 
have the form 

k-q>0 k<0 

= E v/^^w^(i :0+Ev/^^(! -i) (C6) 

k-q>a ^ ^ k<0 ^ ^ 
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will contribute to the spin polarization. By redefining k + q ^ k and k ^ q ^ k in Eq. IC5I and Eq. IC6I separately, we 
obtain the spin polarization on the surface as 

fe>o ^ ^ fc<o ^ ^ 

fc>0 ^ ^ fc<0 ^ ^ 

= E - f^k,t cos{qx)a,_ + ^ Vl - cos(gx)CT^ = ^ 1 - h+ cos{qx)az, (C7) 

fc>0 fc<0 fc 

where /ijl^ is the distribution of the spin-sub band above Dirac point and defined in Eq. IA5I and Eq. IA6I which only 
has finite spin polarization along z direction and will not couple to x direction. 

I 



